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Abstract. We discuss the application of the Mountain Pass Algorithm to the so-called quasi-linear 
Schrodinger equation, which is naturally associated with a class of nonsmooth functionals so that the 
classical algorithm cannot directly be used. A change of variable allows us to deal with the lack of 
regularity. We establish the convergence of a mountain pass algorithm in this setting. Some numerical 
experiments are also performed and lead to some conjectures. 



1. Introduction and results 
The aim of this paper is to find numerical solutions of mountain pass type for the problem 
(1.1) -Au- uAu^ -^Vu= in R^, 

for p > 1 and V > 0. This is the equation of standing waves of the quasi-linear Schrodinger equation 

(i(l)t + + (^A|(^p + 101^" V = in (0, oo) x R^, 
U(0,x) = in R^, 



(1.2) 



where i stands for the imaginary unit and the unknown (j) : (0, oo) x R^ ^ C is a complex valued 
function. Various physically relevant situations are described by this quasi-linear equation. For example, 
it is used in plasma physics and fiuid mechanics, in the theory of Heisenberg ferromagnets and magnons 
and in condensed matter theory, see e.g. the bibliography of [6] and the references therein. The mountain 
pass solutions of the semi-linear equation —Au + Vu = \u\p~^u^ corresponding to ground states for the 
Schrodinger equation 



(1.3) 



\(l)t + + 1^1^" V = in (0, oo) X 
0(0, x) = 0o(^) in R^, 



have been extensively studied in the last decades, both analytically [22] and numerically. On the numerical 
side, the typical tool is the so-called mountain pass algorithm, originally implemented by Choi and 
McKenna [4] (see also [2,3] for a different approach). This works nicely under the assumption that 
the functional associated with the problem is of class on a Hilbert space and it satisfies suitable 
geometrical assumptions. Now, we observe that the functional f : X ^ R naturally but only formally 
associated with (1.1) is 

(1.4) u^S{u):=l [ {1^2u^)\\/ufdx^l [ Vu^ dx ^ / \u\p^^ dx. 

2 JrN 2 Jj^N P + 1 Jj^N 

In fact, taking X as i^^(R^) with > 3, then the functional (1.4) is not even well defined, as it might 
be the case that it assumes the value +oo - (+oo) in the range (TV + 2)/(7V - 2) <p < (37V + 2)/ (TV - 2). 
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In two dimensions, with X = i^^(R^), it is well defined from X to RU {+00} but it is merely lower semi- 
continuous. If, instead, X stands for the set oiu ^ i^^(R^) such that v? G i^^(R^), then it follows that 

(1.4) is well defined from X to R, for every choice 1 < p < (3A^ + 2)/(A^ — 2). From the physical viewpoint 
defining X in this way makes it a natural choice for the energy space. On the other hand, for initial data 
in this space X, currently there is no well-posedness result for problem (1.2) (see the discussion in [6]). 
Furthermore, X is not even a vector space, although (X, d) can be regarded as a complete metric space 
endowed with distance the d{u^v) = \\u — v\\h^ -\- ||Vii^ — \/v'^\\l2 and it turns out that (1.4) is continuous 
over (X, d) suggesting that a possible approach to the study of problem (1.1) could be to exploit the 
(metric) critical point theory developed in [7]. Nevertheless, this continuity property is not enough to 
establish the convergence of a traditional Mountain Pass Algorithm. Actually, it is not even clear what a 
satisfying gradient for £ is. Let us emphasize that we are interested in the convergence of the algorithm 
in infinite dimensional spaces which ensures the convergence of the discretized problem at a rate which 
does not blow-up when the mesh used for approximation becomes finer. Hence, in conclusion, neither the 
mountain pass algorithm can be directly applied for the numerical computation of some solution of (1.1) 
nor, to our knowledge, the current literature (see e.g. [15] and the references therein) contains suitable 
generalization of the mountain pass algorithm to the case of non-smooth functionals, except the case of 
locally Lipschitz functional, which unfortunately are incompatible with the regularity available for our 
framework. On the other hand in [5,6], in order to find the existence and qualitative properties of the 
solutions to (1.1), a change of variable procedure was performed to relate the solutions u G X to (1.1) 
with the solutions v G i^^(R^) to an associated semi-linear problem 

(1.5) -A. = /(.,.), 

^yl + 2r^[v) 

where the function r : R ^ R is the unique solution to the Cauchy problem 

More precisely, u is a. smooth solution (say C^) to (1.1) if and only if v = r~^{u) is a smooth solution to 
(1.5), that is a critical point of the -functional T : i^^(R^) R defined by 

(1.7) v^T{v):=l [ \\/v\^dx ^ / \r{v)\P^^dx^l [ V{x) \r{v)\^ dx. 

2 JrN P + 1 Jj^N 2 J^N 

In addition, as we shall see, the mountain pass values and the least energy values of these functionals 
correspond through the function r. Now, by applying the mountain pass algorithm to T, we can find a 
mountain pass solution v of (1.5). Then u = r{v) will be a mountain pass solution of the original problem 
with a reasonable control on numerical errors. 

We mention that, at least in the case of large constant potentials V and for a general class of quasi- 
linear problems, there are some uniqueness results for the solutions, see e.g. [11]. 

The article is organized as follows. In section 2, we will establish the equivalence between ground 
state solutions for 8 and those for T. Section 3 will recall the mountain pass algorithm and discuss 
assumptions under which its convergence is guaranteed for our problem. In section 4, we will present our 
numerical investigations. We finish by a conclusion giving some conjectures and outlining future work. 

Notation & terminology. We will denote Imf the image of a function f : B i.e., Im/ = {f{x) : x G 
A}. The notation adh(A) stands for the topological closure of the set A. The norm in the Lebesgue spaces 
LP{Q) will be written | -11,^(^7) (^ will be dropped if clear from the context). We will call a (nonlinear) 
projector any function that is idempotent i.e., any function / such that f{x) = x for all x G Im/. 

2. Equivalence between Mountain Pass values 

We noticed in the introduction that, thanks to the change of unknown u = r(v), solutions v G 
H^{R^) n C2(R^) to (1.5) correspond to solutions ueXf) C'^{R^) to (1.1) where X = {u e H^{R^) : 
G i^^(R^)} (see e.g. [5,6]). In this section, we want to show that, in addition, if v is at the Mountain 
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Pass level, then u is at the Mountain Pass level too in a suitable functional setting, see formula (2.3). 
Hence, numerically computing a Mountain Pass solution ^' up to a certain error, yields a Mountain Pass 
solution u = r o V to the original problem (1.1) up to a certain error (involving also the error due to the 
numerical calculation of the solution r to the Cauchy problem in (1.6)). In order to prove this, notice 
that T : H^{R^) defined by (1.7) reads 

T{v) = ]-j \Vv\^dx- j F{x,v)dx 

2 J^N Jj^N 

where we have set F{x^t) := f{x^s)ds with / defined by (1.5). Notice first that, if u = r{v) where 

V e H^{R^) n C2(R^), then ueXn C^iR^). Furthermore, it follows that 

(2.1) S{u) = T{v), 
where E is the action defined by (1.4). Indeed, we have 

S{r{v)) = - [ {l^2r\v)){r')\v)\\/v\^dx ^ / |r(^)|^+Mx+- / V{x)\r{v)fdx 

2 JrN P + 1 JrAT 2 J^N 

= - [ \\/v\^dx — [ \r{v)\P^^dx^- [ V{x)\r{v)\^dx 

2 JrN P + 1 J]^N 2 Jj^N 

= nv), 

thanks to the Cauchy problem (1.6). 

Proposition 2.1. Let v G H^{R^) D C^(R^) be a Mountain Pass solution to problem (1.5)^ that is 

(2.2) T{v) = inf sup r(7(t)), = {l e C{[0,1],H\R^)) : 7(0) = 0, r(7(l)) < O}. 

-feTr te[o,i] 

Then u = r{v) G X D C^(R^) is a Mountain Pass solution to problem (1.1); that is 

(2.3) £{u)= inf sup f(7?(t)), P^ = {7? G C([0, 1], X) : 7(0) = 0, f (7(1)) < O}. 

vers te[o,i] 

Furthermore, u is also a least energy solution to (1.1) for the energy £ (i.e., achieving the infimum of £ 
on non-trivial solutions to (l.l)y). 

Proof It is readily seen that if v e {R^ ) D {R^ ) is a solution to (1.5), then u = r{v) eXnC^iR^) 
is a solution to (1.1). Setting now T-j- := {r o 7 : 7 g Prl^ account of (2.1), we have 

£{u) = T{v) = inf sup T{j{t)) = inf sup £{j{t)). 

On the other hand, we will show that Tj- = Tg, yielding assertion (2.3). In fact, let r] ^Tg and consider 
7 := r-^oT]. Then, r] = ro7 with7 G C{[0,1],H^{R^)), 7(0) = and r(7(l)) = ^(ro7(l)) = £{r]{l)) < 
0. Hence 77 G Pt". Vice versa, if 77 G P7-, there exists 7 G P7- such that 7/ = r o 7. Then, it is readily seen 
that T] G C([0,1],X), 77(0) = r(7(0)) = r(0) = and £{r]{l)) = T{j{l)) < 0. By Theorem 0.2 of [13], 
is a least energy solution. In turn, by Step II in the proof of Theorem 1.3 of [5] the last assertion 
follows. □ 

Remark 2.2. As we will see in the next section, in our case f{v)/v (where / is defined by (1.5)) is 
increasing on (0, +00) and therefore mountain pass solutions v can be characterized by 

T{v) = inf supT{tv) 

instead of (2.2). The numerical algorithm described below finds a local minimum v of v sup^>Q T(tv) 
but there is no absolute guarantee that is a mountain pass solution. Nonetheless, i) is a "saddle point 
of mountain pass type" [15, Definition 1.2] i.e., there exists an open neighborhood V of v such that 

V lies in the closure of two path-connected components of {v ^ V : Tiv) < T{v)}. Since the map 

V r{v) : H^{R^) ^ X is an homeomorphism and £{u) = T(r~^(ii)), the same characterization holds 
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for u := r{v). In conclusion, the above discussion can be thought as an extension of the correspondence 
of proposition 2.1 to saddle point solutions. 

3. Mountain pass algorithms — convergence up to a subsequence 

Let H be a Hilbert space with norm ||-||, E a closed subspace of T : H ^ M a C^-functional and P 
a peak selection for T relative to i.e., a function 

P :n\E ^n-.u^ P{u) 

such that, for any u G 1-L\E^ P{^) is a local maximum point of T on R+i^ ^ E = {tu + e : t > and 
e G E}. The Mountain Pass Algorithm (MPA) uses P to perform a constrained steepest descent search 
in order to find critical points of T. The version we use in this paper is a slightly modified version of the 
MPA introduced by Y. Li and J. Zhou [23] which in turn is based on the pioneer work of Y. S. Choi and 
P. J. McKenna [4]. Let us first recall the version of Y. Li and J. Zhou. 

Algorithm 3.1 (MPA [23]). (1) Choose an initial guess uq G lm{P), a tolerance £ > and let 
n ^ 0; 

(2) if \\VT{un)\\ < £ then stop; 
otherwise, compute 



(3.1) Un+1 ^ . II 

\\\Un - SnVT{Un 

for some Sn G (0, +00) satisfying the Armijo^s type stepsize condition 

T{Un+l) - T{Un) < -^||VT(iin)||||^n+l " ^n||; 

(3) let n ^ n-\-l and go to step 2. 

Conditions under which the sequence {un) generated by the MPA converges, up to a subsequence, to 
a critical point of T on ImP have been studied [24] in the case dim£^ < 00. 

In [21], N. Tacheny and C. Troestler proved that, if an additional metric projection on the cone K of 
non-negative functions is performed at each step of the algorithm, the MPA still converges (at least up to 
a subsequence) and that the limit of the generated sequence is guaranteed to be a non-negative solution. 
In that paper, the authors studied a variant of the above algorithm were the update at each step is given 
by 

(3.2) Un+l = P{PK{Un[Sn])) whcrC Un[s] := - S^^^^^, 

for some Sn G S{un) C (0, +00) with Pk being the metric projector on the cone K. Notice that, with 
this formulation, the projector P only needs to be defined on the cone K. The set S{un) of admissible 
stepsizes is defined as follows: first let 

S'^iun) :={s>0: Un[s] 7^ and T {P Pr {un[s])) - T{un) < -^5||Vr(iin)||} 

and then define S{un) '= S*{un) H (| sup S*{un), +00). Note that the right hand side of the inequality 
to satisfy does not depend on P. Under the following assumptions on the action functional T and the 
peak selection P: 

(PI) P : K \ {0} K\ {0} is weh defined and continuous; 
(P2) inf^^im(p)nK T(^/) > -00; 
(P3) ^ adh(ImPnK); 

(P4) T{Pk{u)) < Tiu) + o{d\si{u,K)) as d\si{u,K) 0; 
they prove that the algorithm generates a Palais- Smale sequence in K. Roughly, using a numerical 
deformation lemma, they show that a step Sn G S{un) exists and that it can be chosen in a "locally 
uniform" way. The trick is to avoid being arbitrarily close to without being "mandated" by the 
functional. Let us remark that the definition of S{un) is natural (but not the only possible one) to 
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force this stepsize not to be "too small". Then, under some additional compactness condition (e.g. the 
Palais-Smale condition or a concentration compactness result), they establish the convergence up to a 
subsequence. If furthermore the solution is isolated in some sense, the convergence of the whole sequence 
is proved. 

Let us also mention that, recently, inspired from the theoretical existence result of A. Szulkin and 
T. Weth [20], the convergence of the sequence generated by the mountain pass algorithm has been proved 
(in some cases, up to a subsequence) even when dim£^ = oo and Q = (see [14]). 

Let us come back to our problem and verify that assumptions (P1)-(P4) hold for the functional T 
given by (1.7) and the projection P defined as follows: P{u) is the unique maximum point of T on the 
half line {tu:t>0}. 

The autonomous case —Au — uAu^ = g{u) was studied by M. Colin and L. Jeanjean [5] who proved 
the existence of a positive solution under the assumptions 

(gO) g is locally Holder continuous on [0, +oo); 
(gl) lim„^o ^ < 0; 

(g2) lim„^+oo|5(w)|/w(3^+2)/(^-2) = when > 3 or, 

2 

for any a > 0, there exists Ca > such that, for all > 0, \g{u)\ < e"^ when N = 2; 
(g3) > such that G(V^) > where G{u) := g{t) dt when TV > 2 (or G(V^) = 0, gi^p) > and 
G{u) < and G (0, V^) when N = 1). 

In our case, this situation corresponds to V ^ constant and g{u) = \u\p~^u — Vu where assumptions 
(g0)-(g3) are verified when p G (3, ^73^)- 

Let us now work with the equivalent problem —Av = r' {v)g{r{v)) where r is defined as in the intro- 
duction and g verifies (g0)-(g3). Let us call f{v) := r'{v)g{r{v)). Recalling that r{v) ^ v when v ^ and 
r{v) ~ ^/v when v +00 [5], (gl) implies that / verifies lim^^o < 0- Moreover, as a consequence 
of (g2), / is subcritical with respect to the critical Sobolev exponent ^^f. Therefore T is well defined 
on H^(W^) and is a local minimum of T. This establishes (P2)-(P3). Because of (g3) there exists a 

such that T{v*) < but, in order to have (PI) and (P4), we need to replace (g3) with the following 
stronger assumptions (see [21]): 

(g4) F{v)/v'^ +00 when v +00 where F{v) := f{s) ds; 
(g5) V f{v)/v is increasing on (0,+oo). 

Let us remark that when the domain is bounded (as in numerical experiments), it is enough to require 
lim^i^o ^ instead of lim^^^o < 0. So we can work with V = in this case. 

For our problem (1.5) with V ^ constant, it is easy to see that (g4) is satisfied. Let us now show 
that property (g5) holds for p > 3. First, let us remark that the derivative of v 
is a positive function times 

(3.3) {k + 2{k - ly) (/ • v) - ((1 + 2r2) • r) . 

For A: = 1, this quantity becomes r^v — (1 + 2r^)r which is negative because r\v) - v < r (see e.g. [5, 
Lemma 2.2]). Thus the map v — Vr(i')/('UY^T+2r^(^) is increasing. As a consequence, it remains 
to prove that f{v)/v is increasing when V = 0. For this, it is enough to show that (3.3) is positive for 
k = p. Using the inequality r\v) - v > (see e.g. [5, Lemma 2.2]), one readily proves the assertion. 

The above arguments imply that the mountain pass algorithm applied to T generates a Palais-Smale 
sequence when 3 < p < ^^^y- Numerically, it is natural to "approximate" the entire space by large 
bounded domains that are symmetric around 0. In the numerical experiments, we will consider 
= {—R/2^R/2)^ . On VLr^ the Palais-Smale condition holds and consequently the MPA converges up 
to a subsequence. This approximation is reasonable for two reasons. First, the solution v{x) on goes 
exponentially fast to when \x\ 00. Second, if we consider a family of solutions {vr) on Qr which is 
bounded and stays away from and we extend ur by on \ Qr^ we will now sketch an argument 
showing that (vr) converges up to a subsequence to a non-trivial solution on R^ (see also [8] where 
authors prove that, for some semilinear elliptic equations —Av + Vv = /(x, v)^ ground state solutions on 
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large domains weakly converge to a solution on R^). The boundedness ensures that, taking if necessary 
a subsequence, vr ^ v. At this point, it may well be that v = 0. Nevertheless, E. Lieb [16] proved that 
if a family of functions {vr) is bounded away from zero in i^^(R^) then there exists at least one family 
(xr) C such that vr{' — xr) ^ 'i;* ^ up to a subsequence. To conclude that v* ^ 0, it suffices for 
example to pick up xr so that /^(^^^ l^R?' 7^ 0- Intuitively, the role of xr is to bring back (some of) 
the mass that vr may loose at infinity. It is then classical to prove that is a non-trivial solution on R^. 
If moreover {vr) is a family of ground states, no mass can be lost at infinity and so vr{- — xr) . In 
this case, v* is a ground state of the problem on R^. Of course, if (xr) is bounded in R^ then vr{- — xr) 
weakly converges to a translation of for any bounded family {xr). In particular, for xr = 0, vr weakly 
converges, up to a subsequence, to a non-trivial solution v. By regularity, vr ^ v in Cj^^(R^) up to a 
subsequence. 

As the moving plane method can be applied to the autonomous case (see [18]), the positive solutions 
Vr are even w.r.t. each hyperplane Xi = and are decreasing in each direction Xi from = to 
Xi = As (vr) is bounded away from 0, the mass of vr is located around 0. Thus vr converges up to 
a subsequence to a non-zero positive solution of (1.5) on R^. Moreover, it is expected that the ground 
state solutions are unique, up to a translation, which would imply the convergence of the whole family 
(vr). This uniqueness result has been proved for V sufficiently large (see [11]). 

The above reasoning can be applied to the sequence generated by the MPA. In Section 4, the numerical 
experiments will provide bounded families {vr) of positive (approximate) solutions on Qr which are 
bounded away from zero (when V 7^ 0). Thus they converge to a solution on R^ as ^ 00. 

For the non-autonomous case {V non-constant), M. Colin and L. Jeanjean [5] proved the existence of 
a positive solution for the equation —Au — uAv? -^V{x)u = g{u) on R^ when 

(VI) there exists Vb > such that V{x) > Vq on R^; 

(V2) lim|^|^+^ V{x) =:Voo eR and V{x) < Ko on R^; 

(gl^) ^(0) = and g is continuous; 

(g20 lim^^o ^ = 0; 

(g30 there exists p < 00 when TV = 1,2 or p < (37V + 2)/{N - 2) if TV > 3 and C > such that 

\g{u)\ < C{1 + \u\P) for any ueR. 

(g4') 3/i > 4 such that < /j.G{u) < g{u)u for any u> where G{u) := Jq g{t) dt. 

These assumptions are clearly satisfied for our model nonlinearity g{u) = \u\p~'^u with 3 < p < ^73^. 
As before, to prove that T possesses the properties (P1)-(P4), we need to replace (g4') with a slightly 
stronger assumption, namely that v /(x, v)/v is increasing on (0, +00) for almost every x G R^ where 
/ is defined by (1.5). Essentially repeating the arguments that we developed for the autonomous case, one 
can show that this assumption is statisfied for our model problem. As a consequence, the MPA generates 
a Palais-Smale sequence. Since the Palais-Smale condition holds for the ground state level (even on R^ 
when V is non-constant), the MPA sequence converges up to a subsequence to a solution of (1.5). 

Concerning the approximation of R^ by large domains, we can argue similarly to the autonomous case. 
Numerically (see Section 4.4), (vr) is bounded away from zero and the peaks of vr are located around 
local minimums ofV. So, the entire mass of vr is not going to infinity. Thus, as in the autonomous case, 
Vr ^ V in Cjq^ with v being a non-trivial solution on R^. 

4. Numerical experiments 

In this section, we compute ground state solutions to problem (1.1) using algorithm 3.1 with the update 
step (3.2) (instead of (3.1)) on problem (1.5). The algorithm relies at each step on the finite element 
method. Let us remark that approximations are saddle points of the functional but, as the algorithm is 
a constrained steepest descent method, it is not guaranteed that they are gound state solutions. Never- 
theless, no non-trivial solutions with lower energy have been found numerically. 

Let us now give more details on the computation of various objects intervening in the procedure. 
As we already motivated above, the numerical algorithm will seek solutions v to (1.5) on a "large" 
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domain := {-R/2,R/2f with zero Dirichlet boundary conditions instead of the whole space R . 
Functions of Hq{Qii) will be approximated using P^-finite elements on a Delaunay triangulation of 
generated by Triangle [19]. The matrix of the quadratic form {vi^V2) ^ Jq^^Vi\/v2 is readily 
evaluated on the finite elements basis. A quadratic integration formula on each triangle is used to 
compute V ^ In^^ \r{v)\P~^^ dx + ^ In^'^i^) \r{v)\'^dx. The function r is approximated using a 

standard adaptive ODE solver. The gradient g := VT{v) is computed in the usual way: the function 
g e Hq{Qr) is the solution of the equation W(p G i^g, {g\(p)fji = dT{v)[(p] or, equivalently, g is the 
solution of the following linear equation 

-Av-f{x,v) innR, 
on dQR, 

where / is defined in (1.5). This equation is solved using the finite element method. 

Since in our case E = {0}, the projector P{u) is easily computed for any u ^ 0: first we pick a 
ti > large enough so that T{tiu) < and then we use Brent's method to obtain the point at which 
t T{tu) achieves its maximum in [0,ti] (alternatively, one could seek t such that dT{tu)[u] = 0). The 
stepsize Sn of (3.2) is determined as follows: we use Brent's method to compute a value of s minimizing 
[0, 1] ^ M : s T{P{PK{un[s]))) ' This choice guarantees that no arbitrary small steps are taken unless 
required by the functional geometry. 

The starting function for the MPA is always (xi,X2) ^ (0.25 - (xi/i?)^) (0.25 - {x2/Rf). The 
program stops when the gradient of the energy functional at the approximation has a norm less than 
10~^. An simple adaptive mesh refinement is performed during the MPA iterations in order to increase 
accuracy while keeping the cost reasonable. The approximate solution is then further improved using a 
few iterations of Newton's method. 

4.1. Zero potential. To start, we study the problem 

(4.1) - Au-uAu^ = \u\P-^u 

in the open bounded domain (—0.5,0.5)^. Running the mountain pass algorithm with p = 4 and = 1, 
we get the results presented in Fig. 1 and Table 1. These experiments suggest that the solutions go to 
as ^ oo, that is when balls become larger. This is not surprising. Indeed, for A/" > 2, the classical 
Pohozaev identity for the semi-linear problem — Av = f{v) on reads 



r r TV - 2 

(4.2) / ^—fiv)v - NFiv) 



dx = 



where F is the primitive of / with F{0) = 0. In our case, / is given by eq. (1.5) with V = and thus 
F{v) = \r{v)\P~^'^ /{p -\- 1). Hence for instance for N = 2^ the condition (4.2) becomes 



/ 



|r(^)|^+^ =0. 



Thus r{y) = or, equivalently, v = {). 

We also repeat the experiments with the larger exponent p = 6 i.e., we consider the problem: 

(4.3) - Au-uAu^ = \ufu. 

The same conclusions can be drawn in this case. 

Then, for p in the range [4, 7], on Fig. 3, we compare the energies and the L'^-norms of the mountain 
pass approximation given previously (our problem with V = 0) with those of the approximation of the 
problem —Au = \u\p~'^u. Equation (4.1) and the latter one are particular cases of the following family. 
Let : R ^ R be the unique solution of the Cauchy problem 
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R = l R = b R=10 R = 30 




Figure 1. Sol. v to (1.5) and u = r{v) to (1.1) on for V = and p = 4. 





R=l 


R = 5 


R = 10 


R = 30 


max V 


415.8 


2.97 


1.17 


0.45 


|Vi;|l2 


875 


5.89 


2.19 


0.80 


maxii 


24.2 


1.77 


0.94 


0.42 


T{v)=S{u) 


78148 


6.75 


1.19 


0.18 


livr(^)ll 


8.2 • 10-^ 


1.4-10-^^ 


1.6-10-^0 


9.6-10-^0 



Table 1. Characteristics of approximate solutions v to (1.5) and u = r{v) to (1.1) on 
Qr for V = and p = 4. 



R = l R = b i?=10 i? = 30 




Figure 2. Sol. v to (1.5) and u = r{v) to (1.1) on Qr for F = and p = 6. 

If 1' is a solution of (1.5) with r = r^, then u = rs{v) is a solution to (1.1) for (5 = 2 and to the Lane-Emden 
equation for 5 = 0. On Figure 3, we have in addition draw the curve for (5 = 1 to shed some light on this 
homotopy. 

4.2. Non-zero potentials. In this section, we study the case V = 10 and p = 4 or 6. We repeat the 
experiments above, see Figures 4 and 5 as well as Tables 3 and 4. As before, we also compare the energies 



MOUNTAIN-PASS ALGORITHM FOR THE QUASI-LINEAR SCHRODINGER EQUATION 9 





R = l 


R = h 


10 


= 30 


max V 


10.10 


2.34 


1.46 


0.80 




18.96 


4.15 


2.52 


1.32 


maxii 


3.59 


1.52 


1.11 


0.70 


nv) = 8{u) 


87.5 


5.00 


1.97 


0.58 


livr(^)ll 


2 • 10-^ 


4.2 • 10-^^ 


3 • 10-^ 


2 • 10-^ 



Table 2. Characteristics of approximate solutions v to (1.5) and u = r{v) to (1.1) on 
VLr for V = and p = 6. 




4567^ 4567^ 



Figure 3. Comparison between (1.1) with y = and the problem —/S.u = \u\p on (—0.5,0.5)^. 

and the L^-norms of the mountain pass approximation given previously (our problem with V = 10) with 
those of the approximation of the problem —Au + lOu = \u\p~'^u for p G [4,7], see Fig. 6. This time, 
when — > oo, the solutions no longer vanish but seem to converge to a non-trivial positive solution on 
R^. This is what is expected in view of the argumentation in section 3. 

We conclude this section by examining the behavior of the solutions as the potential V goes to 0. 
Figure 7 depicts the energy of the solution to (1.1) obtained by the Mountain Pass Algorithm on 
with = 10 for various values of V. It clearly suggests that £{u) ^ as F ^ 0, indicating that the 
ground-state solutions bifurcate from as F approaches the bottom of the spectrum. 



R = l R = b i?=10 i? = 30 




Figure 4. Sol. v to (1.5) and u = r{v) to (1.1) on flji for V = 10 and p = 4. 



4.3. Exponential nonlinearity. Since our numerical experiments are performed in two dimensions, it 
is natural to wonder what shape the solutions are expected to have when the nonlinearity has exponential 
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R = l 


R = b 


i? = 10 


i? = 30 


max V 


417 


11.6 


11.6 


11.6 




877 


20.9 


20.8 


20.8 


maxii 


24.24 


3.87 


3.87 


3.87 


nv) = 8{u) 


79183 


217 


217 


216.7 


livr(^)ll 


2 • 10-^ 


4-10-10 


2 • 10-^ 


1.6-10-^ 



Table 3. Characteristics of approximate solutions v to (1.5) and u = r{v) to (1.1) on 
VLr for V = 10 and p = 4. 



R = l R = h i?=10 i? = 30 




Figure 5. Sol. v to (1.5) and u = r{v) to (1.1) on VLr for ^ = 10 and p = 6. 





i? = 1 


R = b 


10 


i? = 30 


max V 


10.57 


5.98 


5.98 


5.98 


|Vu|l2 


19.56 


9.73 


9.73 


9.72 


maxix 


3.68 


2.68 


2.68 


2.68 


r{v) = 8{u) 


105.1 


47.3 


47.3 


47.3 


livr(^)ll 


7-10-12 


8.8 • 10-^ 


5.6-10-^ 


1.2 • 10-^ 



Table 4. Characteristics of approximate solutions v to (1.5) and u = r{v) to (1.1) on 
VLr for y = 10 and p = 6. 




4567^ 4567^ 



Figure 6. Comparison between (1.1) with F = and the problem —/S.u^l{)u = \u\p 
on (-0.5, 0.5)^ 



growth. To that aim, we consider in this section the equation 

(4.4) - Aii - u/\u^ -^Vu = exp(^/2) - 1. 
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123456789 10 



Figure 7. Energy of the solutions to (1.1) on Qio as V ^ 0. 



The outcome of the numerical computations is presented on Figures 8-9 and Tables 5-6. Unsurprisingly, 
one may see that the qualitative behavior of the solutions is similar to what could be observed before 
with power nonlinearities. 



R = l R = 5 R=10 R = 30 




Figure 8. Sol. to (4.4) on for V = 0. 





R = l 


R = 5 


R = 10 


R = 30 


max V 


7.6 


1.85 


0.31 


0.03 


|V^|l2 


11 


3.76 


0.63 


0.065 


maxii 


3.1 


1.30 


0.30 


0.03 


T{v)=S{u) 


44.4 


2.02 


0.06 


0.0007 


livr(^)ll 


3 • 10-^ 


7-10-^1 


1.7-10-^ 


1.3-10-^ 



Table 5. Characteristics of approximate solutions v and u = r{v) to (4.4) on Qr for V = 0. 



4.4. Non-constant potentials. To conclude this investigation, let us examine the case of a variable 
potential. E. Gloss [12] in dimension N >3 and J. M. do O and U. Severo [9] for = 2 showed that the 
equation 

(4.5) -e'^{Au^uAu'^)^V{x)u = g{u) in 

possesses solutions concentrating around local minima of the potential V when e ^ 0. When the 
potential is constant, this corresponds to R = 1/s and, for a positive potential, the above graphs show 
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R = l R = b R=10 R = 30 




Figure 9. Sol. to (4.4) on for V = 10. 





R=l 


R = 5 


R=10 


R = 30 


max V 


8.04 


6.63 


6.63 


6.62 




10.98 


9.10 


9.10 


9.10 


maxi^ 


3.17 


2.84 


2.84 


2.84 


T{v)=S{u) 


50.60 


41.47 


41.45 


41.40 


livr(^)ll 


3.6-10-'^ 


8.6 • 10-^ 


8.5 • 10-^ 


9 • 10-^ 



Table 6. Characteristics of the approximate solutions v and u = r{v) to (4.4) for V = 10. 



a concentration around the origin on a bounded domain (which suggests that any point will do on R^). 
In this section we have considered (4.5) on (—0.5,0.5)^ with g{u) = \u\^u and the double well potential 

V{x) := 10-8exp(-20|x-c|^) - 5 exp(-30 |x - c |^) 

where c = (—0.2,0.2) and c' = (0.3,-0.2). It is pictured on Fig. 10. Note that the two wells have 
different depths. For e = 0.05, the MPA returns two different solutions (see Fig. 10) depending on the 
initial function. (We have chosen not to display the function v as its shape is similar to the one of u.) 
The one with the lower value of T (hence E) is the solution located around the lower well of V. It is 
obtained by using the MPA with the same initial guess as before. For the right one, the MPA is applied 
starting with a function localized at the other well, namely 

(4.6) (xi,X2) ^ max{0, 0.1 - (xi - 4)^ - (x2 - 4)^}. 

For larger such as e = 0.25, the MPA returns only one solution with a maximum point not too far 
from the point at which V achieves its global minimum, see the rightmost graph on Fig. 10 (the graph 
depicted is the outcome of the MPA with the usual initial function but even using the function (4.6) as 
a starting point gives the same output). For even larger £, such as £ = 1, the solution is very similar to 
the one for R = 1 displayed on Fig. 4. 

5. Conclusion 

In this paper, we tackled the computation of ground state solutions for a class of quasi-linear Schro- 
dinger equations which are naturally associated with a non-smooth functional. A change of variable 
was used to overcome the lack of regularity and a mountain pass algorithm was applied to the resulting 
functional to compute saddle point solutions. 
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Potential V e = 0.05 e = 0.05 e = 0.25 




Figure 10. Potential V and the corresponding sol. u to (4.5) on (—0.5,0.5)^ for p = 4. 



In the autonomous case, we outlined arguments and saw on the above numerical computations that 
the numerical solution ur on converges to a radially symmetric solution u* on the whole space as 
R ^ oo. The existence of was proved by Colin and Jeanjean [5]. The fact that the same numerical 
solution is found with many different initial guesses suggests that the ground state solution on is 
unique, a fact that was proved [11] for V large. The numerical computations also suggest that the set of 
solutions ix* bifurcate from as F ^ 0. 

For the case of variable potential, the numerical computations exhibited the existence of several solu- 
tions of mountain pass type which are local minima of the functional T on the peak selection set Im P. 
The asymptotic profile of these solutions seem to be radial as it is expected from the theoretical results 
on for £ ^ [9, 12]. Interestingly, the multiplicity of positive solutions does not seem to persist when 
€ grows larger. Whether the multiple solutions come through a bifurcation from the branch of ground 
state solutions as e goes to may be the subject of future investigations. 
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